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Abstract. We consider the possibility of experimental verification of vacuum e + e~ pair creation at the focus 
of two counter- propagating optical laser beams with intensities 10 20 -10 22 W/cm 2 , achievable with present- 
day petawatt lasers, and approaching the Schwinger limit: 10 29 W/cm 2 to be reached at ELI. Our approach 
is based on the collisionless kinetic equation for the evolution of the e + and e~ distribution functions 
governed by a non-Markovian source term for pair production. As possible experimental signals of vacuum 
pair production we consider e + e~ annihilation into 7-pairs and the refraction of a high-frequency probe 
laser beam by the produced e + e~ plasma. We discuss the dependence of the dynamical pair production 
process on laser wavelength, with special emphasis on applications in the X-ray domain (X-FEL), as well 
as the prospects for fi + n~ and 7r + 7r~ pair creation at high-intensity lasers. We investigate perspectives for 
using high-intensity lasers as "boosters" of ion beams in the few-GeV per nucleon range, which is relevant, 
e.g., to the exploration of the QCD phase transition in laboratory experiments. 

PACS. 12.20.-m Quantum electrodynamics - 42. 50. Hz Multi-photon processes - 52.38.-r Laser-plasma 
interactions 

1 Introduction A few examples have been discussed of physical situ- 

ations where the Schwinger effect could occur despite the 
Vacuum e+e~ pair creation by a classical electric field is a high critical field strength; e.g., relativistic heavy ion col- 
longstanding prediction in QED [T pp ]. A complete the- lisions [12], neutron stars [T3UH] and focused laser pulses 
oretical description of the effect exists [4j[5l[6j[7] , but there [S] • 

is still no experimental verification. The main obstacle is Since the Schwinger effect is non-perturbative and it 
the high value of the critical electric field strength for pair requires an exact solution of the dynamical equations it 
creation; viz .,E CI = m 2 /e = 1.3 x 10 16 V/cm for electron- is customary to approximate the complicated structure of 
positron casffj]. According to the so-called Schwinger for- a real laser field by a spatially uniform time-dependent 
mula [3], the pair creation rate in a constant electric field, electric field. According to different estimates QHTTJ[l2j 

[T ^ ri ^ fT5 | ri^T J T7T J T8l J T9] the effect of vacuum pair creation 
e 2 E 2 ( irm 2 \ is unlikely to be observable with presently available laser 



S 4^3 exp I | e j5| J ' W parameters. 

However, recent developments in laser technology, in 
is suppressed exponentially when E <C E CT , see Appendix particular the invention of the chirped pulse amplification 
[A] However, a very different situation occurs when the method, have resulted in a huge increase in the light inten- 
field acts only in a finite time interval (dynamical Schwinger sity at the laser focal spot [20l[2Tj . On this basis the Euro- 
effect) [HI5MTU]. In this case, the Schwinger formula, as pean Extreme Light Infrastructure (ELI) project [22] will 
well as its analog for a monochromatic field (Brezin-Itzykson be developed in order to provide radiation beams of femto- 
formula [H]), become inapplicable in the weak field regime, to atto-second duration in the deeply relativistic regime, 

exceeding intensities of 10 25 W/cm 2 . On the other hand, 

Send offprint requests to: David Blaschke construction of X-ray free electron lasers XFEL [2l51P?4"] 

1 We use fi — c — kB = 1 throughout. based on the SASE principle is underway at DESY Ham- 
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burg. Thus an experimental verification of the Schwinger 
effect is coming within reach. 

Under conditions of short duration pulses time-dependent 
effects become important. Therefore in our work |25j we 
use a kinetic equation approach, which allows us to con- 
sider the dynamics of the vacuum pair creation process 
while accounting properly for the initial conditions [81126] . 
Compared to alternative treatments, this approach is es- 
sentially nonperturbative and contains new dynamical as- 
pects, such as longitudinal momentum dependence in the 
distribution functions and non-Markovian character of the 
time evolution. It also takes into account the effects of field 
switching and particle statistics [27lf28l l29] . This approach 
has been applied already to the periodical field case [9J 
with near-critical values of the field strength and X-ray 
frequencies. In particular, it was shown that there is an 
accumulation effect when the intensity of the field is about 
half critical: the average density of pairs grows steadily 
with increasing number of field periods. The method [8] 
also found application in describing the pre-equilibrium 
evolution of a quark-gluon plasma produced in ultrarela- 
tivistic heavy ion collisions at RHIC and LHC [27 | l28 1 l2"9] . 
A characteristic feature of the kinetic approach is the pos- 
sibility of a description of quasiparticle excitations during 
all stages of the external field evolution. 



2 Time-dependent Schwinger mechanism 
2.1 Kinetic approach 

The basic quantity addressed in the kinetic approach [8] is 
the distribution function of quasiparticles in the momen- 
tum representation f(p,t) 4.5,6,7 . The kinetic equation 
for this function can be derived from the Dirac or Klein- 
Gordon equations in an external time-dependent electric 
field by the canonical Bogoliubov transformation method 
[4] or in the oscillator representation approach [26j . This 
procedure is exact but it is valid only for the simplest field 
configurations; e.g., for a spatially-uniform time-dependent 
electric field with fixed direction E(t) = (0,0, E(t)). It 
is assumed that the electric field vanishes at the initial 
time t = to, where real particles are absent ("in- vacuum" 
state). The derivation of the corresponding kinetic equa- 
tion (in the collisionless limit) was presented in Ref. [8J 
and is summarized in the Appendix lAl The result is 



where p is the kinematic momentum and 



at 



dp 



S{p, t) = -A{p, t, t) J dh A{p, h , t) [1 ± 2/(p, t x ) 
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2 J dt 2 e(p,t 2 ,t) 



, (2) 



'■2 

p(t!,t 2 ) =p-e J E(t')dt', 



A(p,t 1 ,t 2 ) = eE(t 1 ] 



e 2 (p,t!,t 2 ) 



> s = 77> 



A(p,t 1 ,t 2 )=eE(t 1 ) g ( Vf v 8 = 0, 
e z {p,ti,t 2 ) 

e(p,t 1 ,t 2 ) = yje 2 ± + p 2 (t 1 ,t 2 ), 



(3) 



wherein m and e are the particle's mass and charge, and 
s is the particle's spin. The created particles are acceler- 
ated by E cxt (t), and the associated currents generate an 
opposing electric field, Ei nt (t). The total field E(t), which 
is finally responsible for particle production, is defined as 
the sum of the external (laser) field E cxt and the self- 
consistent internal field E- mt , which can be found from the 
Maxwell equation 



E int (t) 



e2 2s f dp 



(2n)3j s(p 



{2p„/(p,t) 



2.s 



fdhAfatut) [l±2/(p,i! 



2 / dt 2 e(p,t 2 ,t) 



, (4) 



where e(p) = e(p, t, t) — ^/m 2 + p 2 . 

In general the kinetic equation involves in addition to 
the source term, 5(p,i), its coupling to Maxwell's equa- 
tion, which provides for a field-current feedback typical 
of plasmas, and collision terms. The importance of these 
terms depends on the magnitude of the background field 
and the mass of the produced particles [3D]. For rela- 
tively weak fields we expect the produced-particle number 
density to be small and hence collisions to be rare. We 
therefore neglect the collision term in this work. Quan- 
tum statistics affect the production rate through the term 
[1±2 /] in Eq. ([2]), which ensures that no momentum state 
has more than one spin-up and one spin-down fermion. In 
addition, both this factor and the "cos" term introduce 
non-Markovian character to the system: the first couples 
in the history of the distribution function's time evolution; 
the second, that of the field. 

Equation ^ can be transformed to a system of three 
ordinary differential equations for the functions /, u and 
w 0|, 

df 1 



dt 
du 
~di 
dw 
~dt 



Au, 



A (l±2/)-2e w, 



2e u 



(5) 
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where we have suppressed the arguments. In Eq. ([5|), e = so that for any N > 1 
^/e\ + [<7n — eA(t)] 2 , E(t) = A(t) and the source func- 



tion A is defined in Eq. (|3|) . Another form of this set of 
equations, which includes the Vlasov dynamics, is met, 
e.g., in the treatment of the refraction of a probe laser, 
see Eqs. (|27|) of Subsect. 13.21 The presence of two incom- 
mensurable time scales complicates the problem of direct 
integration of the system: the characteristic time of pair 
creation is \jm and that of field change, 1/v, where v is 
the laser frequency. Furthermore, in our case Eq. ([2]) con- 
tains two small parameters: E <C m 2 and v <C m, but 
we cannot construct any perturbation theory because of 
the memory effects present in the argument of the cosine. 
These memory effects can only be neglected when 



eE m v 
— - < — 

mv 777 



(6) 



but this condition would contradict the quasiclassical con- 
dition for the electric field; i.e., E ^s> v 2 . 

Maxwell's equation, ((4]), can be written in the manner 
of Eq. ©: 



E„ 



J 



~0, 

e 
2^3 

e 
4tt 3 



dp 

£ 



P\\f 



1 

2' 



dp p -l 

£ 



f- 



—u 

2 



s = 0. 



(7) 



In our case, E <C -E'er, the feedback field, E- m t, is negligible 
and / <1 with high accuracy. Under these conditions the 
analytic solution of the kinetic equation is 



f(p,t) 



i J dt 1 A(p,t 1 ,t)exp I 2i J dt 2 £(p, t 2 , t) 



(8) 



The total quasiparticle number and energy density are 
defined as a moments of the distribution function: 



(9) 
(10) 



It is well known j¥j that some of these quantities are diver- 
gent and need regularization, which is discussed in Sub- 
sects. [27J and K2\ 



2.2 Asymptotic expansions and regularization 

In order to regularize divergent expressions for moments 
of the distribution functions, such as (fTO]) . we suggest to 
apply the following procedure [31,32]. Let us expand the 
solutions of Eqs. ([5]) as an asymptotic series in l/e(q), 
where e(q) is given in Eq. J3|): 



f 



51 /n, 

n 

In ~ 1 



W 



E 



00, 



(11) 



N 
n=l 



(12) 



Substituting these expansions in Eqs. ([5]) and separat- 
ing the corresponding orders we obtain for the leading 
terms 



A A (A 

W2= 2-£> U3= 4?' /4= U 



(13) 



Hence the number density integral is convergent but the 
current and energy densities are logarithmically divergent. 
The regularization procedure suggested in [31,32 consists 
in the substraction of the leading terms, Eq. (fT3|) . from the 
integrand: 



jR(t) = 



£ R (t) 



e2 2 



47T 3 

1 

47T 3 



d P P -l 

£ 



f + \(u-u,) 

2 VPu 



2.s 



dpe(f - h 



(14) 



These expressions obey the energy conservation law with- 
out the external field. 

The regularization presented in Eqs. (|14p is equivalent 
to charge regularization in QED [31,32 . Let us write the 
counter-term for the current density in the form 



e 2 E(t) J dp 
\p\>A 



'4e 5 



2.s- 



E(t)U 3 (t), (15) 



where Us(t) is a divergent integral and A is some scaling 
factor with dimension energy. We add the expression in 
Eq. ((TSJ) to both sides of the Maxwell equation and denote 
the regularization constant as Z = (1 + Us)" 1 , then 



Z- l E = -j + V 3 E. 



(16) 



Using the standard definitions of renormalized quantities 
according to the Ward identity: eA = crAr, 

e R = Z^e, E R = Z-^E, (17) 

we obtain the renormalized Maxwell equation 

Er = -3r + E R V 3R . (18) 

Formally the regularization in Eq. (114j) is ambiguous 
because the parameter A is arbitrary. The choice A = 0, 
which was used in |31il32j . is not a preferred value of any 
reason. In our problem we have other energy-scales, e.g. 
(m, eA), which have a definite physical meaning. Choosing 
A ^> to, A 3> eA, is sufficient for the current regulariza- 
tion 

oc 

(eg s ) 2 f dx 

U 3 — R 2 / • 

D7T Z J X 

A 

The inclusion of counter-terms in the region p — > can 
eventually result in some non-physical consequences as, 
e.g., in a negative energy density from Eq. (|14p . 



(19) 
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Fig. 1. Time dependence of the number of produced e + e 
pairs in the volume A 3 for a weak, periodic field, Eq. (|20[l E m = 
3-10~ 5 E CI , A = 795 nm, corresponding to the Thsapphire laser 
|33| (solid line) and for the near-critical field case of an X-ray 
laser [23] with E m = 0.24£ cr , A = 0.15nm (dotted line). The 
values of the residual pair density n r are marked in both cases. 



3 Effects to be tested in experiment 

3.1 Photon production 

In this section we consider the region of field parameters 
already achievable at currently operating laser systems 
|33p34| ; namely, v 2 -C E <C E a , where v is the laser field 
frequency. As criteria for the creation efficiency we use 
the mean density (n) per period and the residual density 
n r , which is taken over an integer number of field periods 
[TU] . Our main result is that optical lasers can generate a 
greater number of pairs than X-ray lasers in the spot vol- 
ume A 3 . This could be observable, e.g., through detection 
of coincident 7 pairs from electron-positron annihilation 
with mean total energy w 1 MeV. 

We consider here a simple model for the laser field that 
can be formed in the focus of two counter-propagating 
laser beams: an harmonic field that acts during z periods 
T = 2717V 



E(t) = E m sin id, 0<t<zT 



(20) 



and vanishes for t < and t > zT. 

Fig. Q] depicts the time dependence of the quasiparti- 
cle pair density, which is generated by the field considered. 
The efficiency of plasma production is shown for the ex- 
amples of a working optical Ti:sapphire laser [33], with 



E,, 



3 ■ 10 Eci, A = 795 nm, and for the planned X- 



FEL at DESY 



with E„ 



0.24EU, A = 0.15 nm. 



The density of e + e~ pairs oscillates with twice the field 
frequency v. The residual density n r , which corresponds 
to an integer number z of field periods, n r = n(zT), is 
negligible in comparison with the mean density (n) for 
optical lasers. The ratio of mean and residual densities is 
well approximated by the relation 



(eE m f 



o 



E 

ZJ 

c 

CD 

E 
o 




(21) 



1 2 3 

Momentum [ p/m ] 

Fig. 2. Momentum distribution (in representative scale) of 
created e + e~ pairs: p stands for transversal momentum (solid 
line) and on longitudinal momentum (dotted line). 



For the case considered [33J (n) ~ 10 7 A -3 and the ratio 
(n)/n r is approximately 10 11 for the field in Eq. ([2"U]) . As 
a consequence, despite the fact that the residual density 
for the X-ray laser exceeds that of the optical laser by a 
large factor, the situation is different regarding the mean 
density: the optical laser produces more pairs in the spot 
volume than the X-ray laser. According to Fig. [T] on av- 
erage there are roughly 10 7 pairs in a volume A 3 , which 
corresponds to a pair density of ~ 10 20 cm" 3 . This esti- 
mate is in agreement with Refs. [351135] . The dense plasma 
exists for the duration of a laser pulse but vanishes almost 
completely after switching off the field. 

This result differs from that obtained with the imagi- 
nary time method [TU] where n(zT) ~ z. That approach 
does not give information about (n). We obtain that the 
mean density of electron-positron pairs is defined only by 
the field amplitude and does not depend on the frequency 
while n r ^ v 2 . Both quantities are proportional to the in- 
tensity of the laser radiation. This results in an accumula- 
tion effect for n r in the near-critical field of an X-FEL [UJ . 
After an integer period number, most of the pairs vanish 
and the residual density becomes negligible in comparison 
with the mean one. The formula in Eq. (|2f [) is inapplicable 
for a pulse shape field, which is a more realistic model of 
a laser beam. In this case the behaviour of n r will depend 
strongly on the parameters describing the pulse shape. 
NB. This is not of concern when solving the equations 
numerically. 

The momentum spectrum of created quasiparticle pairs 
is depicted in Figs. [2] [3J In contrast to the standard as- 
sumption of zero longitudinal momentum for c + -c~ pairs 
[12) , the momentum distribution of the quasipartical pairs 
has a width of the order of m for both transverse and lon- 
gitudinal momenta. The shape of the momentum distri- 
bution is changed drastically when the electric field takes 
a zero value, Fig. [3] a layered structure is formed with 
a characteristic time-scale corresponding to the field fre- 
quency. This peculiar momentum distribution is associ- 
ated with the residual density of e + -e~ pairs. 
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0.000 0.005 0.010 

p[m] 

Fig. 3. Layered structure of momentum distribution (p — |p|) 
at the time of zero field amplitude t — zT/2, associated with 
the residual density n r (in representative scale). 




X [ nm ] 

Fig. 4. Wavelength dependence of the mean density (n) 
(solid line) and the annihilation rate (dotted line) according 
to Eq. J22J for a field strength E = 3 ■ 10~ 5 E CI . 



The quasiparticle plasma created in the weak field case 
can be manifested in various physical effects, such as non- 
linear Thomson scattering [37], damping of electromag- 
netic waves [19] , one and two-photon annihilation [38] , the 
non- linear Breit- Wheeler process [39], etc. As an example, 
we estimate the intensity of two-photon annihilation in the 
plasma volume. The corresponding 7-quanta with mean 
total energy « 1 MeV can be detected outside the focus 
of the counter-propagating laser beams. The production 
rate for this process is 



dN 
dVdt 



J dpidp2cr(pi,p2)/i(pi,t)/ 2 (p2,t) 



x v/(vi - v 2 ) 2 - |vi x v 2 | 2 , (22) 

where v is the particle velocity and a is the cross-section 
for two-photon annihilation [38] 



CT(P1,P 2 ) 



x In 



2m¥(T- 1) 
y/r + Vt - 1 



( r 2 +r _ 1/2 ) 



-(r + l)v/r(r-l) 



(23) 



The t-channel kinematic invariant r is given by 



(pi +P2) 2 
4m 2 



1 



Am 2 



[(£i+£ 2 ) 2 -(pi+p 2 ) 2 ]. (24) 



We have estimated the number of annihilation events 
per laser pulse via Eq. (|2"2"j) with the following parameters: 
pulse intensity I — 10 20 W/cm 2 , pulse duration tl ~ 85 
fs, wavelength A = 795 nm, crossing diameter of laser 
beams « 2.5 /xm [33]. The result is approximately 5 — 10 
annihilation events per laser pulse. The wavelength depen- 
dence of the quantities discussed is shown in Fig. [4] 

We have argued that the simplest model of the laser 
field predicts the creation of a dense quasiparticle plasma 
in the foci of counter-propagating optical laser beams with 
parameters corresponding operating lasers 33,34 . The 



plasma lives during the laser pulse and vanishes almost 
completely after switching off the field. The mean den- 
sity is denned by the field strength and does not depend 
on frequency, reaching values ~ 10 20 cm -3 for a realistic 
field strength of 10 11 V/cm. A manifestation of the plasma 
may be the emission of pairs of 7-quanta with a spectrum 
peaked in the vicinity of a total energy of 1 MeV, with 
an intensity of 5 — 10 events per laser pulse. This would 
be a non-linear transformation of the soft laser photons to 
7-quanta with a frequency ratio of about 10 6 . 



3.2 Refraction experiment 

We consider the possibility of an experimental verifica- 
tion of vacuum pair creation in the focus of two counter- 
propagating optical laser pulses ("photon collider", [231 
140] ) by means of an interference experiment with an ad- 
ditional ultraviolet probe laser beam propagating across 
the plasma region. We suppose that a short-lived electron- 
positron plasma can be created in the focus of the photon 
collider and that the plasma density is sufficiently large in 
order for its Langmuir frequency to become comparable 
with the optical one. The probe laser beam must suffer 
refraction in this nonlinear optical medium, which could 
be detected by an interference refractometer. We perform 
estimates for conditions that apply to the recently con- 
structed petawatt laser Astra-Gemini, with an intensity 
5 x 10 21 W/cm 2 [41]. The effect turns out to be tiny and 
observation would be challenging. 

Experiments using a short wavelength probe laser to 
diagnose the state created in the focus of colliding high- 
intensity laser pulses are discussed extensively in the lit- 
erature [40,42,43 . The polarized QED vacuum created 
in the "photon collider" acts like a birefringent medium, 
with two indices of refraction depending on the polariza- 
tion of the incoming light. This effect, discussed for the 
first time in Ref. [44] , was subsequently investigated in de- 
tail [4311431147] . In Ref. @D] the resulting ellipticity signal 
was calculated, which results from a phase retardation of 
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Fig. 5. Setup for a measurement of the refraction index by the 
probe laser interference method. The short-wavelength probe 
laser beam \ p <C Ao propagates across the direction of the 



high-power beams. 



one of the polarisation directions when a linearly polarized 
x-ray pulse passes through a hot spot generated by high- 
intensity laser beams. It was shown via the Heisenberg- 
Euler Lagrangian [2 that the ellipticity signal is ~ 10~ 
for conditions of the Polaris laser 33J with a probe fre- 
quency 1 keV, which corresponds to the limit of accuracy 
that can presently be obtained with high-contrast X-ray 
polarimeters. This scheme is applied to the case of the 
polarized vacuum generated by two counter-propagating 
beams in Ref . 43J . In that case the probe signal is directed 
across the high-power beam so that, in principle, a refrac- 
tion of the probe beam can be expected. The rotation of 
polarization is estimated to be on the order of 10 -8 rad. 

Unlike schemes designed to measure the polarization 
of the vacuum, we suggest to consider the interaction of a 
probe laser pulse directly with the pulsating e + e~-plasma 
forming a peculiar spatial "diffraction grating" pattern 
or a one-dimensional photon crystal. In the focus of two 
counter-propagating high-power laser beams, the standing 
wave appears for a short time, so that the Schwinger mech- 
anism starts to work close to the wave antinodes whereat 
the plasma "grating" is formed. The refraction of one of 
two split probe beams can be measured by the usual inter- 
ference method, see Fig. The size of the plasma region 
is of the order of the optical wave length and its density 
pulses with twice the optical frequency, so that a stable 
interference pattern can be obtained with a sufficiently 
short probe signal; i.e, t p ~ Arj/2. We use for numerical 
estimates the parameters of the petawatt Astra-Gemini 
laser [41]: A = 800 nm, I = 5 x 10 21 W/cm 2 . 

Detection of the reflected beam would be technically 
easier but the reflection is strongly suppressed in this case 
owing to the absence of a sharp plasma boundary. The 
characteristic scale of inhomogeneity, AL, in the probe 
beam direction is defined by the cross section of the laser 



focus; viz., X p <C AL. This justifies the application of 
geometrical optics when reflection can be neglected [48] . 

Owing to the dependence of the medium's properties 
(e + e _ -plasma) on the electric field of the basic laser, the 
calculation of the refraction index of the probe signal re- 
quires nonlinear optics. However, the mode decomposition 
of the field which is traditional for such problems cannot 
be used in our case because the effect of pair creation is 
a truly nonperturbative one. We therefore consider the 
weak probe field in linear approximation as the basis for 
an exact solution for the high-power laser field. 

To account for the influence of the basic laser's strong 
field on the dielectric properties of the plasma generated 
by it, we use a kinetic equation in the form of Eqs. (fS]). 
We choose the geometry of the probe field as k = (k, 0, 0), 
which is convenient since it is then possible to consider 
the fields of the basic and the probe laser as collinear and 
directed along the axis OZ. All quantities are decomposed 
in a standard way |49] , 

f = fo + 8f, u=u + 5u, v = w + 5w, 
E = E + £E, 
Sf, Su, Sw, <5E ~ exp (— iet + ikr) (25) 

but with an essential difference that the functions fo,uo, wo 
are the exact solutions of the kinetic equation in the strong 
field Eo, instead of equilibrium ones as in plasma theory. 

To describe spatial dispersion we introduce an essential 
assumption that the kinetic equation, Eq. (J2J), is applica- 
ble for the description of small inhomogcncous perturba- 
tions; namely, we add the standard operators v(d/dr) to 
the left hand side of the equations. We now perform the 
substitution 

p = q-eSA, (26) 

i.e., effect a transition to the kinematic momentum related 
to the perturbation field <5E, then 



df 

at 
du 
at 

dw 
~dt' 



v df 
dr 
du 
dr 

dw 

r 

dr 



e<SE 



e<5E 



e<5E 



df 

dp 

du 
dp 

dw 
dp 



1 



Au, 



A(l-2f)~2e w, (27) 



2e u, 



where v = p/eq. 

We substitute here the decomposition in Eq. I|25p and 
compute the first order corrections necessary for a calcu- 
lation of the current. We can then determine the current 
density, Eq. (7}, 



dp 



P]\ S f + 2 £ - L§u 



2tt 3 J s 

and define a permittivity with the help of the relations 

3a = VapEp, 



(28) 



£ a = da/3 H <?al3- 



(29) 
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Substituting the regularized Sf, 8u, we obtain for the 
transverse part of the permittivity 



(tr = 1- 



2e 2 



d 3 p 



n-u.- J : «{<>- - Al) \ (u>- kv) 



dfr 



2e 2 Ur 0p„ 



eqAq dw r 



+e±(lu - kv) 



+ i- 



8st 



a 1 dp 

(w - kv) 2 - Al 



g_L , IdUr 
~2 J r 



2 dp, 



d 3 p 



tt 2 lu J e (n 2 - A, 



A A) 



du r 



-(1-2/,) 



- £o 



dw 



where 



Q 2 = (w - kv) 2 - 4e 2 , 
A = eE ex (t)e±/el 



(30) 



(31) 



Unlike the non-relativistic case [49], the denominator 
(uj — kv) does not generate a pole for all u>, k, but only 
for k/uj > 1. Since our model does not allow an analytic 
calculation of ttr for the whole range of variables, we take 
early advantage of the dispersion relation for transverse 
fluctuations in the oscillations at the stage of calculating 
the momentum integral in Eq., ([30]) : i.e., we exploit the 
condition 



uj 2 e tr 



. 



(32) 



We therefore write k = (w^, 0, 0) and solve the equa- 
tion by iteration with respect to e tr , choosing ej:° = 1 as 
the initial point. Since the pole is absent in this approxi- 
mation, we can integrate over the polar angle, taking into 
account that the optical frequency satisfies w « m and 
hence (fl 2 — Aq) « — 4e 2 for E -C E CT . On the same ba- 
sis it is possible to neglect the second string in the real 
part of Eq. (|30p . Limiting ourselves to the case of a non- 
degenerate plasma: /o -C 1, we have 



= 1 - 



2e 2 



d 3 p 



sx df, 

2 U r 
Zfc Q 



An dw r 



dp, 4e dp„ 



An 2 u 



fc o 
-2 £j 



du r 



e± dw r 

Aq—-^U t - £o- 



dPu 



(33) 



shows that e^} > 



The numerical evaluation of Eq. ([33 
1. Therefore the calculation of the next iteration in Eq. (I30|) 

must be done by respecting a pole at v 2 = l/e^) using the 
known Landau recipe (Dirac identity) [49) . The contribu- 
tion of the J-function to the integral is localized in the 
region of very large momentum p 3> m and can there- 
fore be neglected. The principal value of the integral is 



6x10" 



4x10 
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Fig. 6. The real part of the transverse permittivity calculated 
in second order of iterations (I30I). 
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Fig. 7. The imaginary part of the transverse permittivity 



calculated at leading order of the expansion in the small 
parameters w/m and Aq/eq: 



Re 4? = 1 



(2n) 2 



2e 2 



d 3 p f 2e 

T r I T 2 

fc o I fc o 

AqP^ dw r 
2e 9pu 



A 2 ,-. 2 



2 , 3u r 
Jr +£±- 



3pj_ 
2e 



(34) 



Figs. [51 show that the variation of the refraction in- 
dex of a probe beam in comparison with the vacuum case 
is tiny. It poses a challenge for experimental observation, 
at least for the geometry considered. Other field geome- 
tries need separate consideration owing to the need for use 
of more complex kinetic equations [50l51| . supposing an 
electric field of variable direction. The conclusion obtained 
is in accord with previous estimates [4T)ll43| . 
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4 Extensions towards the Schwinger limit 

4.1 Application in the near-critical case (XFEL) 

A kinetic equation coupled with Maxwell's equation was 
used to estimate the laser power required at an XFEL 
facility to expose the process of QED vacuum decay via 
spontaneous pair production [9j . A 9 TW-peak XFEL laser 
with photon energy 8.3 keV could be sufficient to initi- 
ate particle accumulation and the consequent formation 
of a plasma of spontaneously produced pairs. The evo- 
lution of the particle number in the plasma will exhibit 
non-Markovian aspects of the strong-field pair production 
process and the plasma's internal currents will generate 
an electric field whose interference with that of the laser 
leads to plasma oscillations. 

X-ray free electron laser (XFEL) facilities are planned 
at SLAC [34], namely the Linac Coherent Light Source 
(LCLS) [52], and the XFEL at DESY which emerged from 
the design studies for the e~ e + linear collider project (for- 
merly TESLA now ILC: International Linear Collider) 
[55] , They propose to provide narrow bandwidth, high 
power, short-length laser X-ray pulses, with good spatial 
coherence and tunable energy. It is anticipated that the 
realizable values of these parameters will enable studies of 
completely new fields in X-ray science, with applications 
in atomic and molecular physics, plasma physics and many 
other fields [55] , 

A unique ability of these facilities is to provide very 
high peak power densities. For example, a P = 0.2 TW- 
peak laser at a wavelength of A = 0.4 nm, values which 
are reckoned achievable with current technology [53], can 
conceivably produce a peak electric field strength 




« 
c 

» ,„ " 
■o 10 

0) 

E 

= io " 



E„=0.1E C , X=0.15nm 
E„=0.1 E„ , ^=0.075nm 



E a = VmocP/(ttA 2 ) « 1.2 x 10 16 V/r 



(35) 



Boosting P to 1 TW and reducing A to 0.1 nm, which is 
theoretically possible [51] , would yield an order-of-magnitudc 
increase: E g = 1.1 x 10 17 V/m. Electric fields of this 
strength are sufficient for an experimental verification of 
the spontaneous decay of the QED vacuum 23, 55, 561157], 

A single laser beam cannot produce pairs [16j . (For a 
light-like field F^ U F^ V — and hence the vacuum survival 
probability is equal to one.) Nevertheless, if two or more 
coherent beams are crossed and form a standing wave at 
their intersection, one can hypothetically produce a region 
in which there is a strong electric field but no magnetic 
field. The radius of this spot volume is diffraction limited 
to be larger than the laser beams' wavelength: r CT > A, and 
the interior electric field could be approximately constant 
on length-scales approaching this magnitude. The period 
of the electric field is also determined by A. Hence at an 
XFEL facility one might satisfy the length-scale unifor- 
mity conditions noted above. 

However, the planned laser pulse duration: r co h ~ 100 fs, 
is large compared to the laser period: t 7 ~ 1 as, and thus 
the time-dependence of the electric field in the standing 
wave may materially affect the pattern of observed pair 
production; i.e., vacuum decay might be expressed via 
time-dependent pair production. This possibility can only 



0.0002 0.0004 

timers] 



Fig. 8. Time evolution of the number density for Sets I 
and II, Table [1] In strong fields, particles accumulate, lead- 
ing to the almost complete occupation of available momen- 
tum states. In weak fields, repeated cycles of particle cre- 
ation and annihilation occur in tune with the laser frequency: 
fila = 2.0 x 10 18 27TS- 1 ; i?ib = 4.0 x 10 18 27rs _1 . 



be explored using methods of non-equilibrium quantum 
field theory 11.58,59 or an equivalent quantum kinetic 
theory [8I I271I60] . The latter procedure was employed in 
Rcf. 57], wherein it was shown that pair production occurs 
in cycles that proceed in tunc with the laser frequency. 
While that does not lead to significant particle accumu- 
lation, the peak density of produced pairs is frequency 
independent, with the consequence that several hundred 
pairs could be produced per laser period. 

The proposed XFEL facilities offer a first real chance 
of observing the decay of QED's vacuum, a profound and 
nonperturbative quantum field theoretical effect. Never- 
theless, with a quantum Vlasov equation, one can ask for 
more. This equation yields the time-dependence of the sin- 
gle particle distribution function: f(p,t) := (a£(t) a p (t)) , 
and hence can be used to estimate the laser power required 
to achieve an accumulation of e~ e + pairs via vacuum de- 
cay. Furthermore, in quantum field theory the particle 
production process is necessarily non-local in time; i.e., 
non-Markovian, and dependent on the particles' statis- 
tics. These features are preserved by the source term in 
the Vlasov equation. (The Schwinger source term is re- 



Table 1. Laser field parameters are specified in columns one 
and two: Set I is XFEL-like; Set II is strong. Columns three 
and four describe the density, n maa; (i>), and total number of 
produced particles, iV(t>) = \ 3 n ma x(t>), where i> is the time 
at which the number density reaches its (for weak fields, local) 
maximum. A typical laser pulse length is ~ 80 fs 53 . 





A (nm) 


E (V/m) 




(fm" 3 ) 


N(t>) 


Set la 


0.15 


1.3 x 10 1Y 


4.6 x 10" 


13 


~ 10 3 


Set lb 


0.075 


1.3 x 10 17 


4.6 x 10" 


13 


~ 10 2 


Set Ha 


0.15 


1.3 x 10 18 


7.2 x 10" 


10 


~ 10 6 


Set lib 


0.075 


1.3 x 10 18 


6.4 x 10" 


10 


~ 10 5 
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Fig. 9. Number density calculated with Eq = 0.5 E cv . Solid 
line: solution of Eqs. ([2}, Q; dotted line: solution obtained 
using a low density approximation Eq. ((H). The oscillations 
are tied to the laser frequency, see Eq. (I36I) . 
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Fig. 10. Peak particle number density versus laser field 
strength. There is a striking qualitative change at Eq ~ 
0.25 i?cr, which marks the onset of particle accumulation. 



covered in a carefully controlled weak-field limit [50].) 
Consequently, one can identify the laser parameters nec- 
essary to expose negative energy elements in the parti- 
cle wave packets. The system exhibits this core quantum 
field theoretical feature when the time between produc- 
tion events is commensurate with the electron's Compton 
wavelength. Along with accumulation comes the possibil- 
ity of plasma oscillations, generated by feedback between 
the laser-produced electric field and the field associated 
with the production and motion of the e~ e + pairs; and 
also collisions. (These features are reviewed in Ref. [30j.) 

We model an "ideal experiment," and assume r a = A 
and a nonzero electric field constant throughout this vol- 
ume while the magnetic field vanishes identically. This is 
impossible to achieve in practice and therefore the field 
strengths actually achievable will be weaker than we sup- 
pose. Hence our estimates of the laser parameters will be 
lower bounds. Our model is represented by the set of equa- 
tions Q,® and (HOI)- (Table |T] provides our laser field pa- 
rameters.) There are two control parameters in Eq. ((2|): 
the laser field strength, Eq, and the wavelength, A. We fix 
A = 0.15 nm, which is achievable at the proposed XFELs. 
(NB. By assumption, the volume in which particles are 
produced increases with A 3 whereas the field strength de- 
creases with 1/A: there is merit in optimising A.) Our 
study will expose additional phenomena that become ob- 
servable with increasing Eq. 

The quantitative analysis was performed for two exem- 
plary electric field strengths, the weaker field E = 0.1 E cr , 
which should be obtainable at the proposed XFEL facili- 
ties 34,53 and E — E CT , to provide a strong field compar- 
ison. The results are presented in Table [1] and Fig. [8] The 
dynamics of pair creation is qualitatively changed at field 
strengths ~ 0.5E cr . The detailed investigation of this fea- 
ture is presented in Fig, [9l The results for T = t/X < 100 



are accurately fitted by 

n{T; Eq) = a (Eo) sin 2 2ttT + p(T, Eq) T , (36) 
p(T,E ) = p(E )+p'(E )T. (37) 

One can therefore use p(T, Eq) to quantify the rate of par- 
ticle accumulation. This rate is very small for Eq < 0.2 E CI 
57] and, while noise in n(t) prevents a reliable numer- 
ical determination of p, p 1 , it is nevertheless clear that 
in this case lengthening r co h will not materially increase 
the number of particles produced. It is apparent from 
Fig. [5J however, that the situation is very different for 
Eq = 0.5 E CI . The solid curve is described by Eq. ((36)) . 
with a = 1.2 x 10~ n fm~ 3 , p = 5.4 x 10~ 12 fm~ 3 period -1 , 
p' I p = 0.0033/ period. Clearly, the accumulation rate is 
approximately constant. 

FigureflOlshows the peak particle number density, n(t^), 
where — (4k — 3) A/4 is the time at which E(t) is 
maximal during laser period no. k. The figure displays 
a qualitative change in the rate of particle production 
with increasing laser field strength. (n(k\) ss n(t^) if and 
only if there is significant particle accumulation, otherwise 
n(k\) « n(t>).) 

For Eq < 0.25 E cr there is no significant accumula- 
tion of particles: N(t^ ,E ) w N(t^ , E ), etc., just as ob- 
served in Ref. [57] . However, for Eq > 0.25 E CI there are 
more particle pairs after each successive laser period; e.g., 
Eq = 0.35 E CI brings an order of magnitude increase in N 
over the first 100 laser periods. At such values of Eq one 
is in a domain where numerous electrons and positrons 
produced in a single period are accelerated to relative 
longitudinal momenta that are sufficient to materially in- 
hibit annihilation. This ensures that many pairs remain 
when the next burst of production occurs. Consequently 
N(t^ , Eq) grows considerably with increasing k. Accumu- 
lation means collisions can become important and should 
be included in the kinetic equation if quantitatively reli- 
able results are required. 
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X[nm] 

Fig. 11. Dependence of the mean number of particle- 
antiparticle pairs in the volume A 3 as a function of the wave- 
length for a periodic field with intensity / = 3 • 10 27 W/cm 2 
(E m — 10 15 V/cm) for different particle species: electrons 
(solid), muons (dotted) and pions (dashed). 
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Fig. 12. Time dependence of the number of annihilations in 
the volume A 3 for the same conditions as in Fig. [11] 



4.2 Muon and pion pair creation 

The rapid advancement expected in laser technology, pro- 
gressing toward the Schwinger limit case for electron-po- 
sitron pair creation within the XFEL as well as the ELI 
project [61], stimulates the question about vacuum pro- 
duction of charged particles heavier than electrons; e.g., 
muons and pions. 

The mean pair density per field period is inversely 
proportional to the particle mass, Eq. (j2Tj) . so that the 
mean density of muon and pion pairs is estimated as — 
ro e /207 and — n e /273. In order to achieve for and 
/i^ pairs the same densities as for e pairs, one has to 
increase the laser intensity by 2-3 orders of magnitude. 
According to estimates in Sect. 12.11 for the annihilation 
process of e , a rescaling for pairs will result in a few 
hard (~ 100 MeV) 7-quanta per laser shot. For the de- 
termination of the 27 rate from ^ annihilation we use 
here Eq. ([23]) together with the vacuum cross-section for 
two-photon annihilation of pointlikc scalar mesons [62 



c(Pi,P2) 



(r + iyr 



47rm 2 T 2 ]_ s/t — 1 
(2r-l) 



(r-1) 



In (77 + y/r-1) 



(38) 



where r is given by Eq. ([24]) . 

In Fig. [TTJ we show the dependence on wavelength of 
the mean pair densities and in Fig. [12] the number of two- 
photon annihilations as a function of time for e ± , and 
fi pairs at a field intensity of I = 3 x 10 27 W/cm 2 . 

Here we used (for simplicity of numerical calculation) 
an illustrative model: a harmonic laser field persisting for 
an integer number of periods. As a characteristic quan- 
tity we show the density of created pairs in the node of 
the laser field, which is small but nonvanishing even for 



subcritical field strengths, see Fig. [TT] For a more realistic 
model, such as a laser pulse with Gaussian envelope, it 
can be shown that the number of residual pairs for sub- 
critical field strengths is exponentially suppressed by the 
large factor tot, where r is the pulse duration. Therefore 
the only observable traces of the dynamical Schwinger ef- 
fect in the focus of optical lasers are the products of the 
secondary reactions of short-living quasiparticles as, e.g., 
the two-photon annihilation processes 63,164] or quasipar- 
ticle decays. 



In Fig. [13] we show the rate of 27 annihilation pro- 
cesses per volume for e + e~, it + it~ and [i + n~ pairs as a 
function of the laser intensity. The e + e~ annihilation into 
pairs of 511 keV photons is clearly dominating, but may 
be distinguishable from the other processes because they 
generate hard 7-ray photons of about 100 MeV energy. 
According to this result, the rates from e + e~ annihilation 



at presently available laser intensities of 10 22 W/cm 2 
would correpond to those for annihilation of n + ir' 



fi H pairs at intensities of 10 
within the ELI project [22] . 



26 



m 

and 

10 27 W/cm 2 , planned 



Another method of muon pair creation is considered 
in Refs. [55]. It consists in the coherent collisions of e~ 
and e + produced from a positronium atom by a strong 
laser pulse. This mechanism is an alternative to the con- 
ventional one taking place in non-coherent e~e + colliders. 
In our case we have partly similar conditions when the 
created matter represents a "condensate" of pairs with 
zero momentum [4]. Mathematically this is expressed by 
the condition f(p,t) = /(— p, i), where / and / are the 
momentum distributions of particles and antiparticles. In 
the approximation used this correlated state is conserved 
during the time evolution, which assists the reabsorption 
of all pairs when the electric field disappears. 



D.B. Blaschke et al.: Dynamical Schwinger effect and high-intensity lasers 



11 



10 K 
10 ! 



2 10' 



V 10" 

10' 



10" 




Laser intensity [ W/cm 



Fig. 13. Number of annihilations per unit volume and time 
as a function of the laser intensity for e + e~ (solid), tt + -k~ 
(dashed) and At + /i~ (dotted) pairs. 

5 Heavy ion acceleration in vacuum with 
high-intensity laser beams 

We review here some of our investigations on the possi- 
bility of additional vacuum acceleration of heavy ions in 
laser fields of different configurations [66,67]. The main 
observation is the existence of a big variety of acceler- 
ation modes owing to many fitting parameters, even for 
only one beam of Gaussian shape. For crossed beams there 
are even more options. The relevant variables, such as ini- 
tial velocity or pulse duration, produce an essentially non- 
monotonic dependence of the energy gain that makes the 
search for the most effective mode acceleration very com- 
plex. A threshold level of intensity exists (~ 10 24 W/cm 2 ) 
when the ion moves in the capture mode in one direction. 
The crossed-beams scheme is not less than three times 
more effective than the one-beam scheme in the consid- 
ered range of parameters. However, such a scheme works 
for a definite phasing of beams only, which is difficult to 
obtain at high field intensity. In such a regime other non- 
linear effects (pair creation and vacuum polarization) can 
also be active. 

The idea of particle acceleration by means of laser 
fields was first proposed in Refs. [rJ51f6"9"j and has been in- 
vestigated carefully for a long time, with attention mainly 
to methods using plasma-beam interaction. The laser wake- 
field accelerator concept was introduced in Ref.,[7^ and 
has become widespread in recent years. It is now close to 
practical realization [7T|[72] . A variation of this method is 
the interaction of high- intensity laser pulses with solid tar- 
gets, which has proven to be effective for applications such 
as proton radiography or isotope production and nuclear 
activation [73] . 

The rapid development expected for "table-top" high- 
intensity lasers [6TI174"] opens the possibility for additional 
vacuum acceleration of heavy ion beams and has been 
discussed, e.g., in the context of the upgrade of existing 
accelerators, such as the Dubna nuclotron [75]. In previous 
work [66 we proved that such a scheme should be realiz- 



able at a laser intensity of the order of 10 25 — 10 26 W/cm 2 . 
The presence of plenty of fitting parameters, upon which 
the result depends in an essentially non-monotonic fash- 
ion, makes it complicated to find the most effective accel- 
eration mode. The authors of Ref. [75) have concluded that 
circular polarization is most effective for electron acceler- 
ation. In contrast we have found that linear polarization is 
more preferable for heavy ion acceleration. Furthermore, 
vacuum pair creation, which is initiated at such high field 
intensity [9, 23, 25 , can change considerably the conditions 
for ion beam motion. 

Since a direct calculation of the creation rate for a 
Gaussian shaped field is too complex for the present sta- 
tus of the approach, we perform an estimate with the 
help of the field invariants. The presence of short lived 
plasma matter can cause acceleration of particles in it- 
self even without participation of the laser field, as was 
recently demonstrated in the SLAC experiment [77] . The 
correct description of ion-laser beam interaction at inten- 
sities above 10 22 W/cm 2 must account for non-linear vac- 
uum effects as well as the plasma-beam interaction. 



5.1 One beam scheme 



Following Refs. [76 , 78J we consider heavy ion acceleration 
as the motion of classical point particles in the prescribed 
electromagnetic field 



p = e(E + v x B) 

S 



E 



Ez = ~Tz^^ x + y) s[ll h + ( x ~ y) sin (^ii + < t > p)\ 

[COS0J_ ± COS (<f)± + (f> p )], 



kw 2 
w S 



2w 

s 



(39) 

(40) 
(41) 



B * = 7^2 [( x - v) s[n h ~( x + y) sin (h + &»)]> ( 42 ) 

B X = —Ey, By = E x , (43) 

S = E a cxp { - (2i/r) 4 } exp (- ^J^) , (44) 



where Eq is the field amplitude, (f>u and cf>± are the phases 
of longitudinal and transversal field components: 



arctan 



= 4>q + Tj + arctan 



z x 2 + y 2 



Z C J "' 

.2 



(45) 

i > ( 46 ) 



with rj = ut — kz and z c = kw^/2 being the diffraction 
length, and 



= w 2 [l + (z/z c ) 2 



(47) 



where wo is the minimum spot size of the laser beam at 
focus, 0o is a constant and 4> p is a parameter fixing the 
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Fig. 14. Time dependence of the energy of fully ionized gold 
ions Au 79+ for laser fields of different polarizations. 
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Fig. 16. The energy gain A vs. wavelength A of the laser field. 
For the Phelix laser A = 1053 nm 79 , while for the Astra laser 
A = 750 - 850 nm [HJ. 





10 


c 




o 




01 


8 


o 




Z5 




£= 




> 


6 


CD 




(D 






4 


CO 




CD 




>, 




CO 
1— 


2 


cd 




c= 




LU 








1 2 3 4 5 6 
ty Q [ radian ] 

Fig. 15. Dependence of the ion energy gain A on the initial 
phase, 4>o, for different polarizations: 4> p = (solid line), 4> p — 1 
(dashed line) and 4> p = 0.5 (dotted line). 

polarization of the field (t/> p = 0: linear, <f) p = n/2: circular 
polarization). 

We suppose that at t = to the ion is injected on beam 
axis at a distance zq from the focus with initial velocity vq 
along the axis. The Cauchy problem for Eq. ([55]) . x(io) = 
(0, 0, z ), v(i ) = (0, 0, v ), is solved numerically using the 
standard Runge-Kutta method with basic parameters sets 
specified in Tabled 



Table 2. Basic set of parameters 



I [W/cm 2 ] 


A [nm] 


r [T] 


</>0 


4> P 


w [A] 


5 x 10 25 
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e 


M [GeV] 


v [c] 


to [r] 


zo [A] 




79 


179 


0.968 
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Fig. 17. The energy gain A vs. the intensity of laser field for 
different polarizations. The wavelength A = 1053 nm corre- 
sponds to that of the Phelix laser [79] . 



The main characteristic of interest is the ion's energy 
gain, A, during the period of field action. The choice of ion 
type (fully ionized gold Au 79+ ) corresponds to conditions 
at the Dubna nuclotron [75j where the ion beam has final 
energy 4 GeV/nucleon. 

A picture typical of the ion energy's time dependence 
is depicted in Fig. 1141 The monotonic increase in energy 
indicates the heavy ion capture in the laser field. In these 
conditions the relative phases have a crucial influence on 
the acceleration mode. Figure [15] illustrates the depen- 
dence of the energy gain on the initial phase (f>Q. The most 
significant influence of c/>o is observed for linear polariza- 
tion and the least for circular polarization. The permanent 
character of acceleration is also confirmed also by Fig.[TT51 
where the dependence of ion energy on laser wavelength 



D.B. Blaschke et al.: Dynamical Schwinger effect and high-intensity lasers 



13 




19,0 19,5 20,0 20,5 



x[T] 

Fig. 18. The energy gain A vs. the laser pulse duration at 
small time scale of order of the laser period. 




x[T] 

Fig. 19. Energy gain A vs. the laser pulse duration at large 
time scale; viz., much greater than the laser period. 



is presented. Larger wavelengths entail greater ion energy 
because of the increased interaction length. The depen- 
dence of the energy gain on laser field intensity is similar, 

Fig. mi 

It is interesting that the dependence of the energy gain 
on pulse duration, r, has an essentially non-monotonic 
character both at small time scales, of the order of the field 
period for linear polarization shown in Fig. ll8i and at large 
time scales for circular polarization, shown in Fig. 1191 The 
curve in Fig. ll9l represents simultaneously the envelope for 
rapid energy oscillations with linear polarization. 

The dependence of the energy on the initial ion velocity 
shows similar behavior, Fig. [501 The position of the local 
extrema on this curve can be changed by varying other pa- 
rameters, such as the pulse duration or the field intensity. 
Such a property can in principle be used for decreasing 
the velocity dispersion in the ion beam ("cooling"): maxi- 
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Fig. 20. Energy gain A vs. the initial velocity of the ion. 

mum energy gain is enjoyed by ions within a narrow band 
of particular mean velocities, with little gain at the nodes. 



5.2 Two crossed beams scheme 

The use of special geometrical schemes to increase the 
laser acceleration efficiency has been discussed by many 
authors 80,81,82J. The basic idea is to place the charged 
particles at the crossing point of two laser beams, which 
creates a longitudinally pulling electric field. The result- 
ing field is defined mainly by the phase relations of the 
colliding laser beams. To construct those relations can be 
a difficult technical problem in the case of high-power sys- 
tems [H] . The two laser method is applied successfully in 
plasma accelerators, see also Ref. [53]. 

Let us consider two laser beams propagated according 
to the model of Ref. [37j, Fig.HH 

E x = (E xl + E x2 ) cos 6 + (E zl - E z2 ) sin 0, 
E z = -{E xl - E x2 ) sin 9 + (E zl + E z2 ) cos 9, 
xi.2 = x cos 9 =p z sin 9, 

z\, 2 = ±x sin 9 + z cos 9, (48) 

where Ei(xi,t) and E 2 (x 2 ,i) are defined by Eqs. ([40]) - 
(144[) . To obtain a purely axial field along the z axis it is 
sufficient to set Eqi = —Eq 2 [37] . This provides ID-motion 
of the charge. 

The main characteristic of this scheme is the angular 
dependence of the energy gain, presented in Fig. [22] for a 
basic set of parameters. This picture depends essentially 
on other parameters such as the initial position zq or the 
pulse duration r. For example, the curve A{9) becomes 
smoother by decreasing |zo|. In any case, an optimal value 
of crossing angle exists and is ~ 20°, which correlates well 
with the result of Ref. [82] . The efficiency of the crossed 
angle scheme is about three-times higher than the one 
beam scheme for the range of parameters considered. 
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Fig. 22. Energy gain A vs. the half crossing angle. 



Hence it is natural to suppose that the creation rate is 
negligible for such a field. In contrast to this, the crossed 
beams can produce a field of electric type locally by appro- 
priate geometry and phasing; e.g., in the case of counter- 
propagating beams, when a standing wave is formed. The 
arbitrary crossed angle scheme can only be estimated nu- 
merically, see Fig. [23] We observe that regions with a suf- 
ficiently high value of £ occupy an appreciable part of field 
domain. One can therefore expect the pair creation mech- 
anism to work. The presence of electron-positron pairs 
probably has a positive influence on the efficiency of ion 
acceleration. However, this conjecture must be confirmed 
separately. 



5.3 Background pair creation 

As we have seen, Sect. l3.1[ a homogeneous electric field al- 
ready creates a dense quasiparticle plasma at an intensity 
of the order of 10 22 W/cm 2 . It is therefore necessary to 
account for the possible influence of this effect on the pro- 
cess of heavy ion acceleration. Unfortunately, the direct 
calculation of pair creation in a field of Gaussian shape is 
at present beyond reach. We provide an estimate via the 
field invariants 

E 2 -B 2 , E B, (49) 

supposing that the degree by which these quantities differ 
from zero indicates the probability for pair creation in 
the electromagnetic field. This argument is based on the 
known property pair creation is absent for a plane wave. 

In a sense a Gaussian beam is close to a plane wave 
because the longitudinal component is suppressed by a 



6 Summary 

In this contribution we have described two types of non- 
perturbative QED effects for which the recent develop- 
ment of new, powerful laser technologies has opened up 
the possibility of experimental verification. 

In the first part we reviewed an approach to the de- 
scription of charged particle-antiparticle pair production 
in time-dependent external fields (time-dependent Schwinger 
effect) by a non-Markovian source term in a quantum ki- 
netic equation, which is derived from a quantum-field the- 
oretical formulation, as is summarized in the Appendix. 

As possible experimental signals of vacuum pair pro- 
duction we consider: e + e~ annihilation into 7-pairs; and 
refraction of a high-frequency probe laser beam by the 
e + e _ quasiparticle system, evolving in the focus of two 
counter-propagating optical laser beams with intensities 
from 10 20 -10 22 W/cm 2 , achievable with present-day petawatt 
lasers, and up to the Schwinger limit, 10 29 W/cm 2 , to be 
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reached at ELI. For the highest intensities, the subthresh- 
old production of and 7r + 7r~ pairs should manifest 
itself in observable hard 7 photons from pair annihilation. 

In the second part we summarized explorations of the 
question of whether high-intensity lasers can be used as 
"boosters" for ion beams in the few-GeV per nucleon range. 
We provided estimates for the dependence of the energy 
gain on wavelength, geometry and intensity of the laser 
booster, which show that under optimal conditions a power 
of at least 10 25 — 10 26 W/cm 2 is required in order to 
achieve an energy gain in the range of a few GeV/nucleon. 
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A Derivation of the kinetic equation 

In this section we demonstrate the derivation of a kinetic 
equation with a method of the time-dependent canonical 
transformations [H[H] for particles obeying Fermi statis- 
tics. The analogous derivation for the bosonic case is given 
in [8] We start from the QED Lagrangian 

£ = 0»y (fi^ + ieAJip - m^V - \f^F^ , (51) 

where is the field strength, the metric is taken as 
gf*" — diag(l, —1, —1, —1) and for the 7- matrices we use 
the conventional definition [84] . In the following we con- 
sider the electromagnetic field as classical and quantize 
only the matter field. Then the Dirac equation reads 



(52) 



Using a simple field configuration with the vector potential 
in the Hamiltonian gauge A M = (0, 0, 0, A(t)) and a homo- 
geneous electric field E(i) = (0,0, E{t)), E(t) = -A(t), 
one looks for the solutions of the Eq. (fBTj) in the form 



ij°d + j k p k - e-f 3 A(t) + m 

xx (±) (p,He^, (53) 



where the superscript (±) denotes eigenstates with the 
positive and negative frequencies. Herein the spinors R r 
(r = 1,2) are eigenvectors of the matrix 7° 7 3 satisfying 
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the condition R^R S = 2S rs . The functions ^^(q, t) obey 
the oscillator-type equation 

+ [ w 2 (q J t) + * e A(i)] x (±) (q,<) 







(54) 



where we define the total energy e 2 (q, t) = e]_ + p 2 (t) , 
the transverse energy e 2 — m 2 + q 2 , and the longitudinal 
momentum P\\(t) = q n — eA(t). The solutions X fa)*) 
of Eq. (|54p for positive and negative frequencies are fixed 
by their asymptotic behavior at to = t — > — oo, where 
A(to) = 0. The field operators ip(x) and ^(x) can be 
decomposed by the complete and orthonormalized set of 
spinor functions (l53|) as follows: 



tp(x) 



r.q 



^\x) d± c 



(55) 



The operators 6 qr ,6 qr and <i qr ,<i qr describe the annihila- 
tion and creation of particles and antiparticles and obey 
the standard anticommutation rules. The time evolution 
leads to the mixing of states with positive and negative 
energies and, therefore, non-diagonal terms in the Hamil- 
tonian corresponding to Eq. (|52|) emerge. The diagonaliza- 
tion of the Hamiltonian, which is equivalent to the tran- 
sition to the quasiparticle representation, is performed by 
the time-dependent Bogoliubov transformation 



Mt) = Q q(0 &qr + /3q(t) dt qr , 
dq r (t) = a-q(t) dqr - /9_q(t) b + 



(56) 



qr 



with the imposed condition |a q (t)| 2 + |/3 q (t)| 2 = 1. The 
new operators & qr (i) and d qr (t) describe the processes of 
quasiparticle creation and annihilation. By the virtue of 
Lagrange multipliers, one can find from the equations of 
motion (|54|) that the coefficients in the Bogoliubov trans- 
formation (|56|) are connected via the relations [4] 



d "<" - w e2, * ,qJ "'" ■ 
*<«» - -SfS Q ' (i) e " M<,J " J ' 



(57) 



where the dynamical phase is defined as 

0(q,t o ,t) = / dt'e(q,t') . (58) 

To absorb the dynamical phase it is convenient to intro- 
duce new operators 



Bqr(t) - M*) e- ie ^*<"*) 



(59) 



A,r(*) = dqr(*) e-*^*"'*) 

which obey the anti-commutation relations: 

{i? qr (t),B+ v (t)} = {£> qr (i),£»+ v (t)} = <W S qq , . (60) 

These operators satisfy the Heisenberg-type equations of 
motion 



dt 2e 2 (q,t) 
dD qr (t) _ eE(t)e ± 



dt 



2e 2 (q,t) 



B+(t)+i [H(t), D qr (t)] , 



(61) 



where H (t) is the Hamiltonian of the system of quasipar- 
ticles 

r.q 

- D_ qr (t) D± qr (t)] . (62) 

The first term in the r.h.s. of Eqs. (|6Tj) arises because of the 
transition from the representation (|55j) to the quasiparticle 
one which are unitary non-equivalent. 

Next we consider the evolution of the distribution func- 
tion of particles with the momentum q and spin r defined 
as 

/ r (q,t) = (0 in |6+ (t) & qr (*)|0m) 

= <0 m |B+ (t) B qr (t)|0 m ) . (63) 

According to the charge conservation the distribution func- 
tions for particles and anti-particles are related as f r (q, t) — 
f r {— q,t). Taking differentials in Eq. (f6"3"| with respect to 
time t we have 



df r (q,t) eE(t)sj 



dt 



e 2 (q,t) 



Rc{<P r (q,t)} 



(64) 



Here the function $ r (q, t) = (0i n |D_ qr (t) £ qr (i)|0j n ) de- 
scribes the vacuum production of pairs in the external 
electric field E(t). Applying the equations of motion (|6Tj) . 
one finds 



d$ r (q,t) eE(t)e ± 



dt 



2e 2 (q,t) 



[2/,(q,t)-l] 



2ie(q,i) # r (q,t) . (65) 



The solution of Eq. ([55]) with the initial condition ^> r (q, to) 
may be written in the following integral form 



Mq ,t) = £ -±J d t'^- ) [2fM,t')-i] 

to 

x exp [2i0(q, i', t)] . (66) 

Inserting this result into the r.h.s. of Eq. (|64|) we obtain 
the anticipated kinetic equation 



df r (q,t) _ eE(t)e ± f ^ eE{t')e A 



dt 



2e 2 (q,<) 



£ 2 {q,t>) 



[l-2/ r (q,t')] cos [2 e(q,t',t)\ . (67) 



Since the distribution function does not depend on spin, 
the subscript r can be dropped: f r = /. The substitution 
p = q - eA(t), results in the KE ©. 

The derived equation demonstrates several interesting 
features, such as the dependence on particle longitudinal 
and transverse momentum, the account for spin and statis- 
tics, and the non-markovian character of the time evolu- 
tion. The memory effects are caused by the time integra- 
tion over the statistical factor (1 — 2/) and the non-local 
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cosine function, while the structure of the coefficient A(p) 
defines the shape of the momentum distribution of created 
particles. 

Based on microscopic dynamics, this kinetic equation 
is exact within the approximation of a time-dependent 
homogeneous electric field and the neglect of collisions. 
The source term is characterized by the following features: 

— The kinetic equations |(3J) are of non-Markovian type 
due to the explicit dependence of the source terms on 
the whole pre-history via the statistical factor (1 ± 
2/(p, t)) for fermions or bosons, respectively. The mem- 
ory effect is expected to lead to a modification of parti- 
cle pair creation as compared to the (Markovian) low- 
density limit, where the statistical factor is absent. 

— The presence of the dynamical phases, 8(p,t), under 
the cosine in the integrand (|6"6")) generates high fre- 
quency oscillations. 

— The appearance of such a source term leads to entropy 
production due to pair creation and therefore the time 
reversal symmetry should be violated, but it does not 
result in any monotonic entropy increase (in absence 
of collisions) . 

— The source term and the distribution functions have 
a non-trivial momentum dependence resulting in the 
fact that particles are produced not only at rest as 
assumed in previous studies, e.g. Ref . [85] . 

— In the low density limit and in the simple case of a 
constant electric field we reproduce the pair production 
rate given by Schwinger's formula 

5 cl = lim (2vr)- 3 5 ( d 3 P S(p,t) 

e 2 E 2 ( ttto 2 \ 

= i^ ex H"w)- (68) 

In our approach the electric field is treated as a gen- 
eral time dependent field and hence there is no a priori 
limitation to constant fields. However our result allows to 
explore the influence of any time-dependent electric field 
on the pair creation process. It is important to note that 
in general this time dependence should be given by a self- 
consistent solution of the coupled field equations, namely 
the Dirac (Klein-Gordon) equation and the Maxwell equa- 
tion. This would incorporate back reactions as mentioned 
in the introduction. 

Finally we remark that the source term is characterized 
by two time scales: the memory time 

T"mem T~, (^^) 

eE 

and the production interval 

Tprod = 1/<S>, (70) 

with < S > denoting the time averaged production rate. 
As long as E <C m 2 /e < £j_/e, the particle creation pro- 
cess is Markovian: r mem <C T pro d- 



